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An apparatus is described which can be used for PVT and compressibility measurements on supercritical 
fluids from near room temperature to 600 "C and pressures to 35 MPa. Two separate experimental techniques 
are employed to obtain PVT data over a broad range of the state surface. Burnett expansions are performed to 
generate compressibility factor (or equivalently density) data along a well-behaved supercritical isotherm. A 
series of isochoric measurements is then made to extend the temperature range. Densities assigned to the 
isochores are determined from their intersection with the previously measured Burnett isotherm or gravimetri- 
cally. A computer is used for experimental control and for data logging. Isochoric measurements lasting several 
days can be performed routinely and without operator attention. The apparatus has been tested on propane to 
a temperature of 325 "C. The density data, estimated accurate to ±0. 1 percent, are in excellent agreement with 
other existing data. 
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Introduction 

Perhaps the single most important property of any 
compressible fluid is its PVT behavior. An accurate and 
precise description of the PVT surface can, in principle, 
be used to calculate most equilibrium properties, PVT 
data are also required to analyze experimental data on 
other quantities and to perform engineering design cal- 
culations on fluid handling systems. Wide range PVT 
data for the many technically important fluids at 
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elevated temperatures are, however, often scarce and in 
some cases are nonexistent. 

We describe here an experimental apparatus that can 
be used for PVT determinations on fluids above room 
temperature. Two separate experimental techniques are 
used to obtain data over a broad range of the state sur- 
face (see fig. 1). The well known Burnett [l] 1 method of 
successive gas expansions is used to establish the density 
behavior for the fluid along a well-behaved supercritical 
isotherm. The maximum temperature of such mea- 
surements is, however, limited by such things as the 
compatibility of valve packings and the pressure trans- 
ducer with the elevated temperatures. The measurement 
of pressure versus temperature along a nearly isochoric 
path (pseudo-isochores) in a second experiment is then 
used to extend the temperature range. Densities assigned 
to the pseudo-isochores can be determined from their 
intersection with the previously measured isotherms. 

'Figures in brackets indicate literature references and explanatory 
notes at the end of this paper. 
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Figure 1- Schematic of the P-T 
plane for a typical fluid. Densities 
along the Tb isotherm are ob- 
tained using a Burnett expansion 
technique. P-T measurements 
along pseudo isochores are then 
made to extend the temperature 
range from near critical tem- 
perature la T^. Isochore densi- 
ties are obtained from their inter- 
section with the T B isotherm. 



T B <300°C 



■max 



600°C 



Principal objectives in the design of the apparatus 
were: 1) the apparatus should be capable of mea- 
surements covering a broad temperature and pressure 
range, 2) the determined densities should be indepen- 
dent of the data accuracy for any calibration fluid, and 
3) the apparatus should be automated as much as possi- 
ble to minimize the necessity for operator attention. 

The apparatus is capable of operation from near room 
temperature to about 600 "C at pressures to over 35 
MPa. Densities can be determined by using two com- 
pletely independent techniques, the Burnett method in- 
volving successive fluid expansions, or a gravimetric 
technique involving only the well known density of wa- 
ter used for volume calibration purposes. A small desk- 
top computer is employed for experimental control, 
data logging, and partial real time data analysis. Mea- 
surement runs along isochores taking days to complete 
can be performed automatically and without operator 
attention. 

Experimental 

The experimental apparatuses are illustrated sche- 
matically in figure 2. The Burnett system [shown in a) 
and b) of fig. 2] consists of three spherical cells of differ- 
ent volumes and a diaphram type differential pressure 
transducer with valving as shown. The cells with nomi- 



nal volumes of about 200, 120 and 40 cm 3 allow different 
volume ratios to be selected. These elements are sus- 
pended from a thermally insulated platen which forms 
the top closure to a thermostat ed bath used for tem- 
perature control. This arrangement is similar to that 
described previously in detail [2], 

Mixed liquid bachs are usually preferred and have 
traditionally been used for thermostating in this type of 
experiment. Because of the elevated temperatures of this 
work, however, a commercially available oven with 
vigorous circulation was modified for the purpose. Air 
is normally used as the convective heat transfer fluid 
but, for safety reasons, a continuous gas purge (typically 
nitrogen) is employed for work on flammable samples. 

Even though the heat transfer gas is vigorously 
mixed, substantial temperature differences can exist 
over the oven volume. To minimize these differences, 
the components are encased in a massive aluminum 
block [not shown in a) of fig. 2 for clarity]. The typical 
arrangement is shown in b) of figure 2. Space is pro- 
vided around each component to partially decouple it 
thermally from the block and to permit free flow of the 
heat transfer gas around the components. This arrange- 
ment reduces temperature differences to a few tenths of 
a degree. Further reduction in temperature differences 
is accomplished by using very low power shimming 



376 



Fill-Vac 



PRT 






CZI CP £ZD 




a) 



Fill-Vac 




VSAAAAAAAn mh 



b^M 



WZ, 



^/r/m^v/M,-9/y/>//yA 



i 
mmmmM, 



■/////. ''/////'/,,; // 



b) 




>PRT 



O 



Figure 2 -Schematic of the apparatuses, a) Burnett system showing cells, valving and piping, b) Burnett system showing the arrangement for 
temperature control, c) Isochoric system. In a), V h V lt V % = three spherical cells and DPI=differential pressure transducer. In a), b), c), 
F= forced convection fan and MH = main heaters. In b) and c), SH=shim heaters, PRT = platinum resistance thermometer and 
AB = aluminum block. In c), QBTG = quartz bourdon tube gauge (pressure) and C = condensate cell. 
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heaters located in the air stream above each major com- 
ponent. 

By means of diaphram type differential pressure trans- 
ducers, pressures in the Burnett experiment are mea- 
sured by referencing to oil pressure derived from an oil 
operated commercial pressure balance with a stated ac- 
curacy of 0.015%. 

A spherical configuration was chosen for the cells to 
allow a more accurate calculation of the change in cell 
volumes with pressure. Each cell was constructed with 
an identical ratio of outside to inside diameter to provide 
maximum error cancellation in the calculated volume 
ratios. 

The system employed for the isochoric measurements 
is illustrated in c) of figure 2. This system consists of a 
single spherical cell with nominal volume of 120 cm 3 , 
encased in an aluminum block and suspended from a 
second insulated platen. To change experiments the 
Burnett system is simply replaced in the oven by the 
isochoric system. 

Since the isochoric measurements can be made up to 
the maximum temperature of 600 °C, it is necessary to 
locate the valve and the pressure transducer in a less 
hostile environment. The valve and piping are enclosed 
in a second oven, maintained at the temperature of the 
initial Burnett measurements, where the density has al- 
ready been determined, to allow a more accurate cor- 
rection to be made to the assigned densities along the 
isochores. 

Pressures along the isochores are determined from the 
vibration frequency of a commercial 0-40 MPa quartz 
transducer, calibrated frequently against the primary 
dead weight pressure balance. The transducer is esti- 
mated to be accurate to better than 0.05% full scale. The 
pressure transducer is located outside of the ovens in an 
enclosure thermostated at 40 °C. 

Volumes of the isochoric cell, piping, valve, and pres- 
sure transducer were calibrated by low pressure nitro- 
gen gas expansions from a volume previously calibrated 
by water weighing. This procedure allows isochore 
densities to be determined by weighing the sample after 
a run and provides a completely independent check on 
the densities determined in the gas expansion mea- 
surements. 

Temperatures in both experiments are measured by 
using a platinum resistance thermometer (PRT) cali- 
brated on the IPTS-68 by the National Bureau of Stan- 
dards. The PRT is used to measure the temperature of 
the aluminum block. Temperatures of the components 
are then calculated from the block temperature and 
from any small temperature differences existing be- 
tween the block and the components as determined by 
suitably placed differential thermocouples. 



In use, either apparatus is suspended from a support- 
ing structure. A mechanism for manipulating the ovens 
to enclose the apparatus for an experimental run and to 
retract them to provide access to the components is 
provided. For safety reasons, when working with toxic 
materials, the support structure is enclosed in a manner 
similar to a ventilation hood and is continuously aspi- 
rated by an exhaust fan located on the laboratory roof. 

Computer Control and Data Acquisition 

A computer was used for both experiment control 
and data acquisition. The series of Burnett mea- 
surements necessarily involved manual valve and pres- 
sure balance manipulations. In this case, the computer 
functions primarily as a thermal controller to bring the 
various elements of the apparatus into thermal equi- 
librium at the same temperature. The PRT and thermo- 
couples are interrogated under computer control at a 
programmed time interval, usually about 30-40 s, using 
a scanner of custom design. Scanner relays were pur- 
chased commercially and had a claimed thermal EMF 
of less than 0.02 juV. Potentials are measured with a 
DVM with nanovolt resolution, Primary thermometry 
is performed by taking the ratio of the potential across 
the PRT and the potential across a series connected 
standard resistance of 50 XI, carrying a current of less 
than 0.001 A, using a single DVM for the measurement. 
Two ratio determinations are made for each tem- 
perature measurement with current reversed for each. 
Averaging the two measurements minimizes the errors 
due to residual thermal EMF and tends to cancel DVM 
zero offsets. Temperature differences between the vari- 
ous apparatus components, calculated by the computer 
from the measured differential thermocouple potentials, 
are used to activate the shimming heaters to bring all 
elements to the target temperature. In equilibrium, tem- 
perature differences are typically maintained to better 
than ±0.025 °C at 250 °C. Since corrections are made 
for these temperature differences in the data analysis, 
differences of as much as 0.1 °C are considered toler- 
able. Data are logged onto hard copy and cassette tape 
after equilibrium is achieved, as determined by the oper- 
ator. 

The P-T measurements along pseudo-isochores are 
performed almost completely under computer control. 
Only initial filling of the system and subsequent venting 
or weighing of the sample at the end of a run is required 
of the operator. 

In operation, the computer selects the prepro- 
grammed temperature, sets the oven temperature and 
then determines the block temperature, cell tem- 
perature, and the rate of approach of the cell to the 
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target temperature by interrogating the various thermo- 
couples and the PRT at programmed intervals of about 
30-40 s. When the cell temperature is within 0.1 °C of 
the target temperature, the computer establishes a hold- 
ing mode where the temperature control is continuously 
refined. Temperatures normally oscillate only a few 
hundredths of a degree about the target temperature in 
this mode. After an equilibration time (usually one hour) 
the data are logged onto hard copy and cassette tape, but 
only if programmed limits of temperature deviation and 
temperature and pressure drift rates are satisfied. The 
experiment is automatically shut down after a preset 
maximum pressure or temperature is reached. The sam- 
ple is then manually vented, or it can be condensed into 
a small cell maintained at liquid nitrogen temperature 
for subsequent weighing, 

Isochoric runs on propane at temperatures up to 
325 °C, taking 50-60 h to complete, were carried out 
completely without operator attention. An automatic 
dial-up smart MODEM is interfaced with the control- 
ling computer and allows experimental progress to be 
monitored and controlled by telephone from any remote 
site using a portable terminal or personal computer. Ex- 
periments usually are run continuously including eve- 
nings and weekends. 



Data Analysis 

The usual Burnett experiment involves the expansion 
of the sample gas from an initial volume V x into a second 
evacuated volume V 2 . Pressures before and after the 
expansion are measured. After isolation of V\ and evacu- 
ation of V 2 , the process is repeated successively to low 
pressure. For any particular expansion. 



Pi _ Z t 



■N„ 



where P = pressure, 

Vi+Vz 



z=- 



(1) 



p= density and 



N, 



pRT 

It is assumed here that the temperatures of 



Vi and V 2 are constant and identical and it is recognized 
that the volume ratio N t is slightly pressure dependent. 
A number of ways to analyze such data have been 
suggested [4-6]. One common way is to express the com- 
pressibility factor Z as some function of pressure 



Z=F(P), 



(2) 



Least squares analysis can be applied to the experimental 
data using eq (1) to determine the volume ratio JV at 
zero pressure and the parameters of F (P). If F(P) is a 
virial equation the determined parameters can be identi- 



fied with the virial coefficients provided that the range 
of the data falls within the region where the virial equa- 
tion is valid and the proper number of terms have been 
retained in the virial series. Considerable attention has 
been devoted to the problem of how to determine the 
proper number of terms to retain in the virial expansion 
[4-6], The object of this work, however, is to obtain 
densities. We avoid the problem of the virial coefficients 
and use F(P) only as a fitting function. The method can 
usually be used with data up to a maximum pressure 
corresponding to a density somewhat less than the crit- 
ical density. Densities for higher pressure points can be 
obtained from the lower pressure data using the cell 
constant since 



P/ = (p,+i)iV,. 



(3) 



Even though it is possible to obtain the cell constant 
No from a single fit to the experimental data, it is com- 
mon, and often more accurate, to obtain it from a sepa- 
rate calibration experiment using a highly ideal fluid, 
usually helium. Extrapolation of the pressure ratio 
Pj/Pi + 1 to zero pressure yields the volume ratio directly 
since Z is identically one at zero pressure. A major 
advantage of the Burnett technique is that it is not neces- 
sary to do extensive, individual volume calibrations on 
the system, since only the volume ratio, which can be 
obtained from the expansion data, is required for data 
analysis. 

For this work, the analysis is complicated by the fact 
that small temperature differences are allowed to exist 
between the volume elements composing the system. 
Equation (1) becomes 



Pi+i /=i Pf,i Ti+\ t jZi + i w j 



(4) 



where the system is now considered to consist of n 
volume elements which include the primary volumes V { 
and V 2 plus the volumes of connecting piping, valves 
and pressure transducer all at slightly different tem- 
peratures. The particular volume element is identified 
by the second subscript while the first subscript indi- 
cates the particular expansion as before. In general it is 
no longer sufficient for accurate data analysis to know 
only the cell constant No, the magnitude of all volume 
elements must also be known individually as well as 
certain density derivatives of the state surface. 

For this analysis, the concept of a cell constant is 
retained but now 



JV =S 



iki 



0,1 



(5) 
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where the subscript again indicates P,=0 and all vol- 
ume elements are at the target temperature T. The V itj 
can be expressed in terms of the V 0J as 



Pft, 



V u =Vo,j(l + aAT,j)(\+p J P i ), 



(6) 



where a = volume thermal expansion coefficient, 
/J, = pressure dilation coefficient for volume,/ and the 
AT u = Tij — T are the deviations from the target tem- 
perature T. If the A7*<7*and recognizing the fact that 



H (1 +£,/>) z l+itJ l > 



(7) 



eq (4) can be used to obtain, after some manipulation and 
small approximation, 

AWjm (1 +o A/;-,) Z- #+4-a) 2 ^ AT i+lJ . (8) 

Equation (8) is used to obtain the cell constant N using 
calibration data measured with helium. The first term on 
the right contains measured quantities. The second term 
is of order 10 -4 A^ which allows the use of calculated 
values for the individual volume ratios (Voj/Vo.0 with- 
out introducing significant error in the value of JVo ob- 
tained. 

Using eq (6) in eq (4) the expression describing any 
expansion becomes 

Pi Z «Tja % tW l+j8AH + qAr f+u \ m 

P l+I l+aAT i , 1 +p l P i Av A Z l+l ,jT l+UJ J'™ 

where the very small terms in a/J have been ommitted 
for clarity but which are included in the actual data 
reduction programs. For AT^T it is assumed that 

dZ- ■ dZ- 

— =f* =~ and the compressibility factor Z,- at the target 

oJ/j a J 

temperature T is related to Z t j by 



i i J_32,| AT 

1+ 2, ar|, Ar ' J 



(10) 



Equation (9) can then be written, again after some ma- 
nipulation, as 



Pi 



i+t 



Zj + \ 



1+ ±jZi 
^Z, dT 



ma 



A+P i+i B+aC 
+0 l P i +aAT i 






where 



A=N + 



J_^£i: 



2 ^d AT, 



pM-l 9 r \rJ=l Vo,t 



(+!./> 



(11) 



(12) 



and 



£=2 ^0j-t 



Y^ 



1 apy-H i 



ar 



iJ^Ar, +!i , 



^T^Ar,. +u+ - 



1 9p/+i 



aJ 



2^ 
w-i "o.i 



A?W 



(13) 



(14) 



For Ar^O this reduces to the isothermal Burnett situ- 
ation. The Aq in eq (12) is the largest term and is deter- 
mined from the helium calibration runs. The remaining 
terms in (12), (13) and (14) are of order 10~ 4 No or 
smaller and are treated as small correction terms for 
which calculated values of the (Vaj/V 0i i) are again used 
and for which it is sufficient to use only approximate 
values for the derivatives. Taking the logarithm of (11) 
to get 



laP, - lnP i+ , = InZ, - lnZ /+ 1 + lnJV 4- 
[corrections < 10~ 4 ] 



(15) 



yields the working equation which was used, together 
with eq (2) in the analysis of the lower pressure Burnett 
data. Densities at the higher pressures in the Burnett 
measurements are obtained from the lower pressure 
densities of that run using the analogue of eq (3), which 
for this apparatus becomes 

p,=p t+ iN' lt (16) 

where A 7 ",' is an effective volume ratio given by 

1 apy+i 

\.,f, v., " J '+ij 

(17) 



\T i „Jl_^C£±A I y ZJ±Li AT 



pi 3T 



Ar u 



Data for Propane 

Numerous PVT data for propane up to 325 °C have 
been published. To verify the performance of the appa- 
ratus and to confirm the accuracy of the data analysis, 
limited measurements on propane to 325 °C have been 
made and have been compared with the literature val- 
ues. Propane samples were research grade certified to be 
of 99.99% purity. 

The two largest cell volumes were used for the expan- 
sions. Four expansion runs at 250 °C using helium were 
first made from a maximum pressure of about 5 MPa to 
a minimum pressure of about 0.13 MPa in order to de- 
rive the cell constant N . A value of N = 1 .60702+ .0001 
resulted from the analysis of pressure ratios using eq (8) 
and this value was used in the propane analysis. 

Three expansions runs on propane at 250 °C, consis- 
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ting of a total of 30 pressure observations, were made 
[7], Pressures ranged from a maximum of about 35 MPa, 
corresponding to a density of about 1.6 times the critical 
density, to a minimum of about 0.14 MPa. 

The data are separated into two groups for analysis. 
At lower pressures, corresponding to densities less than 
the critical density, the data from all runs could be ana- 
lyzed simultaneously using an expression for the com- 
pressibility Z in terms of pressure [8]. Since it is the 
nature of Burnett experiments to concentrate the data at 
the lower pressures, only the highest pressure point on 
each run could not be handled in this manner. 

For data below 17 MPa the compressibility factor Z ■ 
was represented as 



Table 2. Densities along the 250 "C isotherm obtained from the 
Burnett measurements and gravimetrically (the latter marked with*). 



= 1+2 5, 
i= l 



<(0 



(18) 



where P = 17 MPa and was used with eq (15) for Fitting 
purposes. Data were fitted, in the sense of least squares, 
by minimization of the difference between the left and 
right sides of eq (15). For pressures below 1.3 MPa, 
differences were weighted by P~ x to reflect the de- 
creasing relative accuracy of the pressure balance at the 
lower pressures [8], Other data were assigned a weight 
of 1. However, the results of the fit were insensitive to 
the weighting. For simplicity, derivatives appearing in 
the small correction terms of eq (15) were evaluated 
from a recently formulated correlation of propane data 
[9]. In the absence of derivative data, however, an iter- 
ative procedure can be used with equally good results. 
In this case a best estimate of the derivatives, based on a 
corresponding states approach, for example, could be 
used, or the correction terms could be ignored in first 
approximation. At the completion of the experiments, 
when the isochoric data are available, the derivatives 
can be estimated and successive iterations performed. 
Coefficients for eq (18) are shown in table 1. Burnett 
density results are tabulated in table 2 along with the 
densities obtained gravimetrically during the isochoric 
measurements. 

It must be emphasized that eq (18) is simply a fitting 
function, and its coefficients cannot be identified with 
the virial coefficients, although the coefficient of the 

Table I. Coefficients of equation (18). 



B, 


= -4.14084758401E- 


-01 


B 2 


= 




B, 


= 




B 4 


= 4.01931812102E- 


-01 


B 5 


= -2.20024107867E- 


-01 


B 6 


= 9.14412474575E- 


-03 


B, 


= 4.22392090629E- 


-02 



P (MPa) 


p (mol/dm 1 ) 


P (MPa) 


p (mol/dm 3 ) 


34.614 


7.857 


3.712 


0.937 


33.953 


7.784* 


3.501 


0.879 


28.960 


7.241* 


2.962 


0.734 


25.118 


6.714* 


2.390 


0.583 


22.311 


6.252 


2.250 


0.547 


20.322 


5.864 


1.894 


0.456 


20.277 


5.843* 


1.521 


0.363 


16.261 


4.891 


1.430 


0.341 


12.922 


3.891 


1.200 


0.284 


12.719 


3.827* 


0.960 


0.226 


12.185 


3.649 


0.902 


0.212 


10.384 


3.044 


0.755 


0.177 


8.536 


2.421 


0.603 


0.141 


8.083 


2.271 


0.566 


0.132 


6.922 


1.895 


0.473 


0.110 


5.674 


1.507 


0.354 


0.082 


5.364 


1.413 


0.221 


0.051 


4.566 


1.179 







linear term is in reasonable agreement with published 
second virials for propane at 250 °C [10,11], 

Densities for the highest pressures on each of the 
three runs were calculated from the next lower density 
of that run using eq (16). 

Pressure measurements were made as a function of 
temperature at fixed total amounts of fluid using the 
isochoric system to extend the temperature range. Pres- 
sures were measured along 10 pseudo-isochores cov- 
ering a density range of 1.7 to 7.7 mol/dm 3 . Obser- 
vations were made at 10 °C intervals up to 150 °C and at 
25 °C intervals up to 325 °C and to a maximum pressure 
of about 34 MPa. For isochores intersecting the 250 °C 
isotherm below 8.5 MPa the densities assigned to the 
pseudo-isochores were calculated from the intersection 
pressures using eq (18). 

Since the data are concentrated at the lower pres- 
sures, the density data above 8.5 MPa were too widely 
spaced to permit accurate interpolation. To fill in these 
gaps several isochoric runs were terminated by con- 
densing the samples into a cell held at liquid nitrogen 
temperature for the purpose of weighing. Sample 
weights were determined to better than ±0.01 gram and 
densities at 250 °C were calculated using the volume 
calibration data for the isochoric system. Uncertainty in 
densities determined gravimetrically is estimated at bet- 
ter than ±0.1% due primarily to the uncertainty in vol- 
ume calibrations. These densities are tabulated in table 2 
along with the Burnett densities. Densities for pressures 
above 8.5 MPa, determined by both the Burnett and 
gravimetric methods, were used in a polynominal fit of 
pressure versus density, which in general represented 
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the densities to better than 0.05%. This function was 
used to calculate the densities assigned to the pseudo- 
isochores intersecting the 250 °C isotherm above 
8.5 MPa. 

In calculating the densities along the pseudo- 
isochores, corrections were made for fluid residing in 
the piping, valve, and the pressure transducer as well as 
for the effects of thermal expansion and pressure dilation 
of the cell. Data along the pseudo-isochores are 
presented in table 3. 

Density data along several isotherms are compared 
with the literature data [11-14] in figure 3. Density devi- 
ations are plotted relative to an equation of state for 
propane [9] used as a baseline. This was done to elimi- 
nate the need for multiple interpolations when com- 
paring data at slightly different temperatures. Absolute 
deviations in figure 3 reflect only the ability of the equa- 
tion of state to reproduce the experimental data. Excel- 
lent agreement between data sets is indicated, however, 
by the fact that all the experimental data lie in a band of 
about ±0. 1% which is well within the combined experi- 
mental errors. The excellent agreement between closely 
spaced densities along the 250 °C isotherm obtained in 
the Burnett measurements, the gravimetric mea- 
surements and from the literature data verify the per- 
formance of the experimental apparatus and confirm an 
estimated accuracy of about ±0.1% for the Burnett 
data. 
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Table 3. Propane densities along pseudo-isochores. 



P (MPa) 


T(°C) 


p (mol/dm 3 ) 
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7TQ p 


(mol/dm 1 ) 


5.191 


99.999 


7.314 


5.606 


119.994 


3.855 


6.707 


109.988 


7.308 


6.175 


129.995 


3.853 


8.260 


119.996 


7.302 


6.741 


140.000 


3.850 


9.831 


130.004 


7.299 


7.296 


149.998 


3.848 


11.416 


140.000 


7.294 


8.678 


174.998 


3.843 


13.003 


149.938 


7.288 


10.041 


199.994 


3.836 


17.003 


174.979 


7.275 


11.388 


225.003 


3.830 


21.015 


199.993 


7.262 


12.719 


249.987 


3.824 


24.993 


224.992 


7.250 


14,040 


275.034 


3.817 


28.960 


249.993 


7.239 


15.341 


299.905 


3.811 


5.570 


1 10.003 


5.899 


16.635 


324.943 


3.805 


6.587 


119.998 


5.894 


7.470 


129.999 


5.734 


7.620 


129.999 


5.890 


8.467 


140.023 


5.730 


8.662 


139.992 


5.888 


9.468 


150.001 


5.728 


9.712 


150.016 


5.884 


11.983 


175.014 


5.718 


12.345 


175.000 


5.873 


14.498 


200.018 


5.708 


14.981 


199.999 


5.863 


17.009 


225.033 


5.698 


17.613 


224.996 


5.853 


19.504 


249.980 


5.689 


20.227 


249.985 


5.843 


21.985 


275.002 


5.680 


4.180 


89.991 


7.875 


24.449 


300.021 


5.671 


5.950 


100.002 


7.868 


26.891 


325.023 


5.662 


7.768 


109.998 


7.861 


4.866 


119.988 


2.653 


9.614 


1 19.994 


7.858 


5.217 


129.989 


2.651 


11.477 


130.004 


7.852 


5.564 


140.030 


2.649 


13.350 


140.007 


7.845 


5.905 


149.996 


2.648 


15.230 


149.991 


7.839 


6.744 


175.023 


2.644 


19.941 


175,006 


7.825 


7.566 


200.039 


2.640 


24.646 


200.015 


7.812 


8.374 


225.034 


2.636 


29.319 


225.001 


7.799 


9.169 


249.968 


2.632 


33.953 


250.005 


7.787 


9.955 


274.960 


2.628 


4.804 


100,000 


6.783 


10.734 


299.980 


2.624 


6.082 


110.002 


6.778 


11.504 


325.041 


2.619 


7.397 


119.996 


6.772 


5.134 


119.995 


3.013 


8.735 


129.999 


6.770 


5.549 


130.011 


3.011 


10.087 


139.996 


6.765 


5.957 


140.011 


3.009 


11.449 


150.021 


6.760 


6.359 


149.976 


3.007 


14.869 


174.997 


6.747 


7.354 


174.989 


3.003 


18.298 


199.996 


6.735 


8.323 


199.998 


2.998 


21.721 


225.006 


6.724 


9.282 


225.023 


2.994 


25.118 


250.006 


6.713 


11.160 


274.984 


2.984 


28.498 


275.021 


6.703 


12.084 


300.032 


2.980 


31.853 


300.032 


6.693 


12.996 


324.977 


2.975 


6.020 


120.012 


4.805 


3.722 


110.000 


1.779 


6.779 


130.006 


4.802 


3.935 


120.001 


1.778 


7.539 


140.011 


4.798 


4.145 


130.004 


1.777 


8.298 


150.001 


4.795 


4.352 


140.009 


1.776 


10.192 


175.012 


4.789 


4.556 


150.001 


1.775 


12.073 


199.970 


4.781 


5.058 


174.999 


1.772 


13.946 


224.994 


4.772 


5.550 


199.996 


1.770 


15.801 


250.005 


4.764 


6.035 


224.995 


1.767 


17.642 


274.943 


4.756 


6.512 


249.994 


1.765 


19.462 


299.904 


4.749 


6.984 


275.029 


1.762 


21.282 


324.917 


4.741 


7.449 


299.971 


1.760 


4.443 


99.998 


3.860 


7.911 


325.020 


1.757 


5.031 


109.999 


3.858 
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Figure 3— Comparison of density 
data from various sources: 
this work, D 

Beattie (ref. 12), O 
Thomas {ref. 1 1), & 
Warony (ref. 13), V 
Tiechman (ref. 14) O 
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